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Abstract. We suggest a new signature for dark matter annihilation in the halo: high 
:nergy positron line radiation. Because the cosmic ray positron spectrum falls rapidly 
vith energy, e + ’s from halo WIMP annihilations cam be a significant, clean signal for 
rery massive WIMP’s (> 30 GeV). In the case that the e + e - annihilation channel has an 
ippreciable branch, the e + signal should be above background in a future detector, such 
is have been proposed for ASTROMAG , and of potential importance as a dark matter 
signature. A significant e + e - branching ratio can occur for neutralinos or Dirac neutrinos. 
3igh-energy, continuum positron radiation may also be an important signature for massive 
leutr&lino annihilations, especially near or above the threshold of the W + W~ and Z° Z° 
innihilation channels. 
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Most of the matter in the Universe is dark, 1 and if f l > 0.15, that dark matter must 
be non-baryonic, as primordial nucleosynthesis constrains Qb to be less than 0.15. 2 There 
axe strong theoretical arguments (e.g., inflation, structure formation, and the temporal 
Copemican principle) favoring fi = 1, in which case non-baryonic matter must account 
for 90% or more of the material in the Universe. There are numerous, well motivated relic 
particle candidates for the dark matter; 1 among them, the neutralino, the axion, a light 
neutrino species, or a heavy neutrino species. The nature of the ubiquitous dark matter is 
certainly one of the most important questions facing both particle physics and cosmology. 

The intriguing hypothesis that relic WIMP’s comprise the dark matter is being tested 
by a number of different and complementary experimental approaches. 3 There are acceler- 
ator searches for supersymmetric partners to the known particles, v mass and oscillation 
experiments, double beta experiments (which double as WIMP ionization detectors), and 
direct searches for the relic particles themselves (axions, magnetic monopoles, particle cold 
dark matter). In addition, there are efforts to search for the annihilation products of mas- 
sive relics (weakly-interacting, massive particles, or WIMP’s) that accumulate in the sun 
and earth 4 (high energy i/P’s) and for the annihilation products of WIMP’s which reside in 
the halo 5 (p’s, 7’s, and e + ’s, including 7-line radiation 5 ). In this paper we will discuss high 
energy (> 30 GeV) positron line radiation, and to a lesser extent high-energy continuum 
positron radiation. 5 Because the cosmic ray positron spectrum falls so rapidly with energy 
above ~ 10 GeV (roughly as E~ 3i , see below; for comparison, the 7-ray spectrum falls 
only as E~ 2A , or so), this is a particularly interesting and potentially promising signature 
for high mass WIMP’s. While there are many uncertainties underlying both the astro- 
physics and the particle physics associated with the problem, the signature, especially a 
positron line, appears promising enough to pursue, as it could be seen with future detectors 
proposed for the ASTROMAG facility. 7 

Cold, thermal particle relics are particle species that were once in thermal equilibrium 
and whose relic abundance arises because their annihilations froze out when the tempera- 
ture of the Universe was ~ 1/20 of their mass. 8 For such relics their abundance today is 
related to their annihilation cross section by 

fl ~ 10 -35 cm 2 / < <r|v| >ann 

where < <r\v\ >ann is the thermally-averaged annihilation cross section (evaluated at 
T ~ m/20). As is now well- appreciated, the present abundance of such a relic is inversely 
proportional to the annihilation cross section. For a value of order 10 -38 cm 2 , the relic 
particles provide closure density. 
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There is every reason to believe that cold particle dark matter will find its way into 
the halos of spiral galaxies, including our own, when they form. It is believed that the halo 
density in our galaxy has a roughly constant value, about 0.3 GeV cm -3 , from the center 
of the galaxy out to about 10 kpc (~ the core radius of the halo), and well beyond this 
decreases as the radial distance squared.® The halo material inferred in spiral galaxies from 
flat rotation curves, in some cases measured out to 3 times the distance where the light has 
has all but disappeared, contributes at least IIhalo — 0.03. 1 Since there is no convincing 
evidence for a rotation curve which turns over, the total mass contained in spiral galaxy 
halos has yet to be determined, and could be as great as J"1halo = 1* We will assume for 
now that the local WIMP density is that of the halo, so that the local number density of 
WIMP’s is 

n ~ 10 -2 rrijQ cm -3 

where 77130 = m/30 GeV. Later we will consider the possibility that WIMP’s only contribute 
a fraction of the halo density. 

High-energy positrons created by WIMP annihilations will accumulate in the halo for 
a time of order 10 7 yrs, the estimated containment time for cosmic ray e ± ’s in the galaxy, 
before they diffuse out of the galaxy. 10 In addition, as they propagate they slowly lose 
energy, the dominant losses at the energies of interest being synchrontron radiation and 
inverse Compton scattering off the 2.75 K background radiation and background starlight. 
The effect of energy loss will be discussed below. While there are uncertainties in the 
estimate of the containment time, it is certain to be greater than the light travel time across 
the halo (only ~ 30,000 yrs). Thus the flux of positrons builds up over a containment 
time, and the integrated line flux is given by 

* _ — — — ( c T g)BR (!) 

Here r is the containment time for positrons in the halo, / is the ratio of the actual WIMP 
mass density to the fiducial value of the halo density we have assumed (0.3 GeV cm -3 ), 
BR is the branching ratio to the e + e~ annihilation channel, and g is a geometrical factor 
which we will discuss shortly. [Note the accumulation effect enhances the positron flux by 
a factor of order r/30,000 yrs; there is no similar enhancement for 7 rays produced by 
WIMP annihilations.] For further discussion, it is convenient to re-express the expected 
positron line flux relative to some canonical parameters, as follows: 
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where Tr = rj 10 7 yrs, BR-i = BRj 10 -1 , and < <r|v| >ann= 3 < <r|u| >_ 36 
xlO -2 ® cm 3 s” 1 (in cgs units). 

The cosmic ray positron flux has been measured up to energies of almost 30 GeV. For 
energies of order 1-10 GeV the positron flux is ~ 10% that of the electron flux; for energies 
of 10-30 GeV the fraction rises monotonically to ~ 25%. 11 The cosmic ray electron flux 
itself has been measured up to energies of about 2 TeV. For energies greater than about 
10 GeV the differential electron flux is given (to within a factor of 2) by 11 

&T-1&E- ~ 0.07(£_/GeV)- 3 ' 3 cm" 2 sr -1 s -1 GeV" 1 (3) 

If as a crude estimate and guide we take the positron flux to be about 5% of the electron 
flux, this then implies an extrapolated differential positron flux of 

dT+ / dE+ ~ 4 x 10 -8 (.E+/30 GeV) -3,3 cm -2 sr -1 s -1 GeV" 1 (4) 

[While there is a paucity of 7-ray data above energies of a few GeV, for comparison, the 
extrapolated, diffuse 7-ray flux is roughly 10 -9 (£7^/30 GeV) -2 * 4 cm -2 sr -1 s -1 GeV -1 .] 

The conventional explanation for origin of the positron flux is that it results from 
interaction of primary cosmic rays (p’s and 4 He nuclei) with nuclei in the ISM producing 
Jf-mesons, 7r-mesons, and p - mesons whose decays produce positrons. 12 The theoretical 
expectation, which to be sure depends upon a purely theoretical model for e + production 
and propagation, is an e + /e~ flux ratio of ~ 10% at energies of 0.3-1 GeV, decreasing 
above energies of a few GeV to a value of 3-5%. 11 In light of the ideas discussed here, it 
is interesting to note that the measured e + /e~ flux ratio actually rises from the expected 
value of ~ 10% at an energy of 1 GeV to ~ 25% at the highest energies measured, a value 
that is ~ 5 times that predicted from conventional sources (see Fig. 1). 

Since the extrapolated positron flux falls as E~ 3 i while the predicted flux from WIMP 
annihilations falls only as m -2 , the prospects for its detection become better with in- 
creasing energy (assuming for the moment that the accumulation time r is energy inde- 
pendent). The positron line from WIMP annihilations is expected to be very narrow, 
A E + ~ mv ~ 0.03m 3 o GeV. However, because of the energy resolution of proposed detec- 
tors (few %) 7 and line broadening due to energy loss (see below), we have expressed the 
extrapolated differential positron flux per GeV. For the canonical values used as normal- 
izations above, the positron line radiation from WIMP annihilations starts to dominate the 
extrapolated positron flux at an energy of ~ 20 GeV. Now let us discuss the astrophysical 
parameters r, g, and / in more detail. 

Due to synchrontron and inverse Compton energy losses the WIMP-produced positron 
line (which is of negligible intrinsic width) will be broadened. The energy loss of a cosmic 
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ray e + (or c ) is given by 13 
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dE 4 E 2 , 

dt ~ 3mi <TT ^ + Pm ^ 


( 5 ) 


where <tt = 0.665 x 10 -24 cm 2 is the Thomson cross section, ~ 7r 2 T 4 /15 = 0.27 eV cm -3 
is the energy density in the 2.75 K background (we neglect the subdominant and position 
dependent starlight contribution), = B 2 /8ir ~ 0.22 eV cm ~ 3 B 2 ( B = 3 B 3 x 10~ 6 G), 
and we have taken the rma average of the sine of the angle between the magnetic field 
and momentum to be i/2/3. The first term accounts for energy loss due to inverse 
Compton scattering and the second to synchrontron radiation. Eqn. (5) can be written in 
a more suggestive form, 

dE E 2 

(6a) 


dt 

tae(Eq) = 


tae(Eo)E 0 

1.2 Gyr 


(£ 0 /GeV)(l + 0.81B2) ^ 

Here t&e(E o) corresponds to the energy loss timescale for an e* of energy Eq, t^e(Eq) = 
—Eo/(dE/di). The field strength of the galactic magnetic field is ~ 3 x 10"® G; moreover, 
our halo population of WIMP-produced positrons may be exposed to a smaller rma field 
strength (however, little information exists about the magnetic field of the galaxy outside 
the disk). 14 Noting these uncertainties we will adopt 2?a = 1. 

We will use the following very simple, spatially-homogeneous model to estimate the 
broadening effect of energy loss on the positron line: a delta function source of positrons 
at energy m and strength a = n 2 < <r\v\ >ann (gf 2 )ER/ 4ir; a diffusion time r ~ 10 T yrs, 
and energy loss given by, —dEjdt = E 2 /rAB(m)m. [Note we are assuming that WIMP- 
produced positrons are not accelerated by any processes in the ISM, and thus can only 
lose energy.] Based upon this simple model the partial differential equation governing the 
differential energy flux is 13 


8 dF+ . . * — - r 

= aS(m) — + 


1 dF+ d( E 2 dT+ \ 
r dE 8E \mr^£(m) dE ) 


at dE 

It is simple to find the following analytical solution for the steady state flux: 

d?+ n 2 < <t\v\ >ann ( gf 2 )(cr)BR mr r _„ _ , mi 

lE = 4w H 2 ex P [ r (l — m /E)\ 


r — 


_ TAg(m) 


(7) 


(8a) 


( 86 ) 


r — (m/30 GeV)pr 

valid for E <m. For energies E > m, the flux of course vanishes. Note that the first term 
in the expression for the predicted flux is just the previous expression for the integrated 
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line flux, cf., Eqn. (1), and that the shape of the line is controlled by the ratio r = energy 
loss time / containment time. 15 The predicted differential flux rises up to an energy m, 
and then sharply drops to zero; the width of the broadened line is, A E+/m ~ r -1 . While 
energy losses do significantly broaden the positron line, the expected sharp drop off for 
energies E > m is a very distinctive signature. Moreover, the predicted e + /e~ flux ratio, 

d^ + /d£ 3.5% mj£(gf 2 ) < <r\v\ >_ 36 Bi2_i(E/30GeV) 13 exp[r(l - m/E )] 

has an even sharper shape (see Fig. 1). 

Next, we mention a possible geometrical concentration effect. The population of 
WIMP-generated cosmic ray positrons is produced throughout the halo of the galaxy. 
In the usual models of cosmic ray propagation, cosmic rays are magnetically confined to 
some region in the galaxy of comparable size to that of the disk (or the disk plus bulge), 
and escape by diffusion, or other processes such as a galactic wind. 12 We note that the 
expected scale height of conventional cosmic rays, a few 100 pc, is very small compared 
to the scale of the halo, about 10 kpc. 12 If cosmic rays are most efficiently trapped within 
the disk (because of larger magnetic field strength and shorter diffusion lengths), then the 
density of halo WIMP annihilation produced positrons may be significantly higher in the 
disk. (This, of course, is where the detectors are). We have accounted for this possibility 
by the geometrical factor g\ on naive geometrical grounds one might expect g to be of 
order the ratio of the volume of the halo to that of the disk, or about 30. 

Finally, it is also possible that the local halo density is greater than the canonical value 
which we have chosen, perhaps by a factor of 2. On the other hand, if IIwimp is less than 
about 0.03, it may be that WIMP’s contribute only part of the halo density, the simplest 
expectation being a fraction of order flwiMp/0.03, implying that / 2 ~ (^wimp A^ -3 )' 
Using the fact that < <r|v| >_ 3 #~ ^wiMP> i* f°U°ws that / 2 < <7-|vj >_ 3fl ~ 10 3 flwiMP 
for Q wimp ^ 0.03. We see that the maximum value of / 2 < cr\v\ >_ 38 ~ 30 obtains 
for ttwiMP — &HALO — 0.03, where / could be order 1 and < a\v\ >_ 3 s order 30. Of 
course, for values of fiwiMP less than ~ 0.03, / 2 < <t\v\ >_ 36 decreases as fiwiMP* On an 
editorial note, we wish to remind the reader that relic WIMP’s from the early Universe 
are interesting and important even if they do not contribute closure density or even the 
halo density. Recall that the fraction of critical density contributed by the microwave 
background is only about 10 -5 , and the microwave background radiation is certainly very 
interesting! 

Clearly there are substantial, irreducible astrophysical uncertainties in our expectations 
for the positron flux. It seems fair to say, however, that values for gf 2 Vj as large as 10 3 
are conceivable. 
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Up to this point we have not been explicit about the identity of the relic WIMP. We 
have done so in part because the positron line signature is generic to any cold particle 
relic. Now we will consider some specific expamples; first, the Dirac neutrino. For a relic 
Dirac neutrino species, 0 ~ (3GeV/m) 2 , so that a 10 GeV Dirac neutrino species only 
contributes about 1% of closure density. 1 ® In addition, the branching ratio to the e + e“ 
channel should be of order 3%. Using our previous formula for a species which contributes 
less than halo density we find that 

/ 2 < *\o\ >—38— lOmfo 2 
T+ ~ 2 x 10 -fl -^-cm -2 sr -1 s _1 

m 30 

Note that because of the dependence of the cross section on the mass, for Dirac neutrinos 
the flux decreases as m ~ 4 (at least as long as m <■ M\y)- Thus the positron line loses 
relative to the extrapolated positron spectrum as one goes to higher energies. In addition, 
the results of the UCSB-Berkeley-LBL Germanium double beta decay experiment seem to 
exclude a Dirac neutrino in the mass range of 15-1500 GeV. 17 We use the phrase u seem to 
exclude,” as the authors of this paper have assumed that the Dirac neutrinos contribute 
the full halo density, irrespective of their mass. Here we have taken the more realistic 
view that they contribute only a fraction of it; whether or not their results exclude this 
possibility is not clear. 

A more promising and well-motivated cold particle relic is the neutralino. 18 For the 
simplest particle physics models of the neutralino, the e + e~ annihilation channel is severely 
suppressed, with a branching ratio of order 10~ s , or so, making it of interest only if the 
astrophysical parameters are favorable. The reason for the suppression is easily understood 
as follows. 19 The neutralino is a self-conjugate (Major ana) fermion. In the s-w&ve annihi- 
lation channel the spatial part of the incoming neutralino wave function must symmetric; 
and so, to insure that the overall wavefunction is antisymmetric, the spin part must be 
antisymmetric. Thus the incoming state has zero angular momentum. For most standard 
processes that contribute to neutralino annihilation chirality is conserved, so that massless 
fermions and antifermions come with opposite handedness. Therefore the angular momen- 
tum along the axis of the outgoing fermions is one, which precludes the s-wave. Thus 
for massless outgoing fermions the s-wave amplitude is zero; for massive fermions it is 
proportional to ntf/m. P-wave annihilation is not so suppressed, but it is proportional to 
the relative velocity of the incoming neutralinos, squared. Generically then, the neutralino 
partial annihilation cross section to the fermion-antifermion channel is given by 

< <*\v\ >ANN~ G 2 P m 2 [(mf/m) 2 + v 2 ] 


(9a) 

(96) 
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Note that in the early Universe, when neutralinos decouple, v 2 ~ 6 T/m ~ 1/3 is not 
so terribly small, and p-wave annihilation is not badly suppressed. Thus unfortunately, 
the cross section that determines the relic density is unsuppressed. 20 In the halo, where 
v 2 ~ 10 -5 , p-wave annihilation is suppressed and neutralino annihilations proceed mainly 
through the heaviest fermion, usually the bottom quark. In this case the branching ratio 
relevant for our estimate for the positron line is order 10 -5 . There are, however, more 
cormlicated supersymmetric models where the right- and left-handed selectron masses 
are 'qual. In these models, if there is mixing between the right- and left-handed 
self. there is an additional contribution to neutralino annihilation into e + e“ which 

is proportional to the mass splitting, but which should otherwise be unsuppressed. 21 In this 
case the e + e~ branching ratio can be very substantial, making the positron line signature 
extremely interesting. 

While we have focused on two specific examples of WIMP’s for which the positron 
line radiation could possibly be interesting, we wish to emphasize the generality of our 
suggested dark matter signature. Aside from astrophysical uncertainties and the mass of 

the WIMP, our estimate only depends upon the e + e - branching ratio (and Owimp in 
the case that Qwimp is not unity). Unless that branching ratio is suppressed for special 
reasons (as can be the case with the neutralino), one expects the branching ratio for a 
generic WIMP to be of order 10%, making the positron line signature of very general 
interest. 

Finally, we should comment on the continuum positron radiation which arises from 
neutralino annihilations in general. Such radiation has been mentioned by other authors; 5 
here we wish to emphasize its possible importance for heavier WIMP’s, and especially the 
existence of additional annihilation channels which can be important for heavier neutrali- 
nos, namely the W + W~ and Z°Z° channels. Continuum positron radiation arises from 
the decays of the neutralino annihilation products. For orientation, consider first the usual 
channel, 6$. The 6 and b decay to charmed quarks and a virtual W\ the decay of the 
virtual W + produces a positron with a branching ratio of ~ 8%; the average energy of the 
positron is ~ 1/4 the neutralino mass. [Of course, the virtual W can also produce r’s and 
fi's which ultimately produce positrons, albeit of degraded energy.] Secondary decays of 
the charmed quarks, and of their decay products, can produce additional positrons. The 
hardest positrons will be those produced in the initial b quark decay, and as we argued 
above it is the most energetic positrons that are of the greatest interest. 

When the neutralino mass exceeds half that of the W (~83 GeV), a new channel for 
neutralino annihil ation opens up: a real and virtual W. Further, when the neutralino 
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mass exceeds that of the W, the channel to 2 real W s also opens up. These channels are 
not suppressed by any consideration of chirality. They can easily compete with or even 
dominate the annihilation into fermions, particularly if winos are lighter than sleptons and 
squarks. About 8% of the time the W decay will produce a positron and neutrino. Around 
threshold the W momentum is not very large so that the decay is like that of a W at rest, 
producing a positron which carries away of order 40 GeV. In addition, the positrons are 
produced in a 2, rather than 3, body decay. Thus the fraction of energy carried away by 
the positron should be larger than from b decay, and the spectrum should be sharper. [A 
similar discussion applies to the Z°Z° channel, although the branching ratio to e + e~ is 
only 3%.] Therefore, the positron signal from neutralino annihilations in the mass range 
near the W mass or above may too produce a detectable positron signal. 

The identification of the composition of the ubiquitous dark matter in the Universe is 
a most important question facing both cosmology and particle physics. Quite correctly, a 
wide range of experimental approaches are being pursued. In this Letter we have pointed 
out one more potential dark matter signal: positron line radiation produced from WIMP 
annihilations within our own halo. While the existence of a detectable signal is by no means 
assured, if found it would be decisive evidence for dark matter of a definite mass. This is 
an additional motivation for experiments to measure the high-energy electron and positron 
cosmic ray spectra, and to be sure that such experiments have adequate resolution. 


Note Added: D. Seckel has called to our attention an unpublished manuscript by A.J. Tylka 
and D. Eichler 23 which discusses in detail the spectrum of cosmic ray positrons expected 
from halo photino annihilations, and also mentions the possibility of positron line radiation 
from Dirac neutrino annihilations in the halo. 


It is a pleasure to thank D. Mueller, M. Kamionkowski, G. Giudice, K. Griest, and 
D. Seckel for valuable conversations and comments. This work was supported in part by 
the DoE (at Chicago). 
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Figure Captions 

FIGURE 1 — The predicted positron fraction, e + /(e + + e“), and the existing experimental 
data (from Ref. 11). In calculating the predicted positron fraction we have taken m = 25 
GeV, r = 1, and gf 2 = 20, values which make the comparison to the existing data intrigu- 
ing! In addition, we have assumed a contribution to the positron fraction from conventional 
sources of the form, e + /(e + + e“) = 0.02 + 0.10(E/GeV) -0-5 , which is consistent with the 
models discussed in Ref. 12. 
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